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Direct detection of Weakly Interacting Massive Particle (WIMP) candidates of Dark Matter (DM) is 
studied within the context of a self-consistent truncated isothermal model of the finite-size dark halo 
of the Galaxy. The halo model, based on the "King model" of the phase space distribution function 
of coUisionless DM particles, takes into account the modifications of the phase-space structure of the 
halo due to the gravitational influence of the observed visible matter in a self-consistent manner. The 
parameters of the halo model are determined by a flt to a recently determined circular rotation curve 
of the Galaxy that extends up to ~ 60 kpc. Unlike in the Standard Halo Model (SHM) customarily 
used in the analysis of the results of WIMP direct detection experiments, the velocity distribution 
of the WIMPs in our model is non-Maxwellian with a cut-off at a maximum velocity that is self- 
consistently determined by the model itself. For our halo model that provides the best fit to the 
rotation curve data, the 90% C.L. upper limit on the WIMP-nucleon spin-independent cross section 
from the recent results of the CDMS-II experiment, for example, is ~ 5.3 x 10^* pb at a WIMP mass 
of ~ 71 GeV. We also find, using the original 2-bin annual modulation amplitude data on the nuclear 
recoil event rate seen in the DAM A experiment, that there exists a range of small WIMP masses, 
typically ~ 2 - 16 GeV, within which DAMA collaboration's claimed annual modulation signal 
purportedly due to WIMPs is compatible with the null results of other experiments. These results, 
based as they are on a self-consistent model of the dark matter halo of the Galaxy, strengthen the 
possibility of low-mass ( ^ 10 GeV) WIMPs as a candidate for dark matter as indicated by several 
earlier studies performed within the context of the SHM. A more rigorous analysis using DAMA 
bins over smaller intervals should be able to better constrain the "DAMA regions" in the WIMP 
parameter space within the context of our model. 
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I. INTRODUCTION 



One of the most favored candidates for the Dark Matter (DM), that constitutes roughly a quarter of the mass-energy 
density of the Universe [l|, is some kind of Weakly Interacting Massive Particles (WIMPs) with masses typically 
between a few GeV and a few TeV. Such particles are predicted in many models of physics beyond the Standard Model 
of particle physics such as models involving supersymmetry or higher (spatial) dimensional theories 043 • Owing to 
their weak interactions, such particles would "freeze out" in the early Universe at a characteristic temperature scale 
of ^ MeV, and would be present in the Universe today as 'cold' dark matter with a calculable abundance that can 
naturally explain the measured average cosmological density of dark matter in the Universe. Also, the measured 
rotation curves of spiral galaxies, such as that of our own, the Milky Way, out to and beyond their visible edges, can 
be explained naturally under the hypothesis that the visible galaxies are embedded in much larger, roughly spherical, 
dark matter "halos" . The true extents and the masses of such halos are, however, unknown. 

Following early suggestions 043 (sse @ for a review) , several experiments worldwide are currently engaged in attempts 
to directly detect these hypothetical WIMP DM particles comprising the DM Halo of the Milky Way (the Galaxy, 
hereafter) in terrestrial laboratories. These direct detection experiments attempt to detect nuclear recoils resulting 
from interaction of WIMPs with a suitably chosen target material in underground laboratory detectors with low 
backgrounds. 

The DAMA/Nal [i| and DAM A/LIBRA experiments (hereafter together simply referred to as the "DAM A" 
collaboration) have reported an annual modulation signal in their event rate at a confidence level of 8.2a based 
on their combined data They interpret this as evidence for WIMPs, ascribing the modulation to the periodic 
variation of the flux of WIMPs passing throu gh the detector on earth caused by the orbital motion of the earth 
around the sun. The CDMS-II collaboration has also reported two candidate events in their "signal region" 
with an estimated probability of ^ 23% of observing 2 or more background events in that re gion , for the estimated 
level of the background. Several other experiments such as CRESST [l^, XENONIO [Hill, COUPP [H and 
PICASSO [3, for example, have reported only null results so far. These results provide important constraints on the 
WIMP mass and WIMP-nucleon interaction cross section. The question of compatibility of the DAMA results [l^l 
with the negative results from other experiments has also been the subject of several recent studies which include the 
earlier overlooked effect of "ion channeling" (see below) on the analysis of the DAMA results • 



In order to derive constraints on (or to determine) the mass and interaction cross section of the unknown WIMP 
DM particles from the results of these direct detection experiments, one requires the density and velocity distribution 
of the WIMPs in the solar neighborhood as crucial input parameters. These parameters are a priori unknown. In 
the customary analysis the DM halo of the Galaxy is assumed to be described by a single component isothermal 
sphere [20| with a Maxwellian velocity distribution in the Galactic rest frame given by 
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Here w = |v|, p(x) is the density at the location x, and (w^)^/^ is the velocity dispersion. The justification for this 
Maxwellian form comes from the work of Lynden-Bell [2l| who argued that evolution of systems of collisionless particles 
(such as WIMPs) under the process of gravitational collapse is governed by the process of "violent relaxation" whereby 
the observationally relevant 'coarse-grained' phase space distribution function (DF) of the system rapidly relaxes, 
due to collective effects, to a quasi-Maxwellian stationary state which is a steady-state solution of the collisionless 
Boltzmann equation.^ 



The DM density in the solar neighborhood is usually taken to be Pdm.© ^ 0.3 ± 0.1 GeV/ cm^, following the 
analysis of the observed velocity distribution of stars transverse to the Galactic disk near the sun, as sug- 
gested by Oort and extended by Bahcall A recent "model-independent" analytical study Q quotes 
Pdm,o = 0.43 ± 0.11 ± 0.10 GeV/ cm'^. Also, a recent attempt to directly extract this density from the results of 
cosmological large N-body simulations yields a value of Pdm,© ~ 0.37 GeV cm~^ [2^. The last Reference also finds a 
DM particle velocity distribution that can in general be called "quasi-Maxwellian" . albeit with significant deviation 



The 'fine-grained' DF may, however, never reach equihbrium. 



3 



from the pure Maxwellian form especially at the high velocity end where it drops off more sharply compared to the 
pure Maxwellian. 

The value of the velocity dispersion, (u^)-'^/^, the single parameter that characterizes the Maxwellian velocity distri- 
bution (dJ of the DM particles, is usually determined from the relation [lO], {v"^)^^^ = -^/fwcoo, between the velocity 
dispersion of the particles constituting a single-component self-gravitating isothermal sphere and the asymptotic value 
of the circular rotation speed, Vc^oo, of a test particle in the gravitational field of the isothermal sphere. Neglecting 
the effect of the visible matter on the DM halo, and assuming Vc.oo ~ Vc.q ~ 220kms~^, where Vc.q is the measured 
value of the circular rotation velocity of the Galaxy in the solar neighborhood, one gets {v^)^^^iso ~ 270kms~^. ^ 
This isothermal sphere model of the DM halo with a value of the DM velocity dispersion (w^)d{i — 270kms~^ and 
local value of the DM density Pdm,© — 0.3 GeVcm"^ is what is often referred to as the Standard Halo Model (SHM). 
This SHM is generally taken as a sort of benchmark model of the DM halo of the Galaxy, and the results of various 
WIMP DM search experiments are mostly analyzed within the context of this SHM. 

Whereas the SHM would suffice for the initial design of the experiments, it is necessary to consider improved models 
to analyze the results of the experiments, which are becoming progressively more sophisticated. The reasons for 
this are many: (1) It is well-known that the isothermal sphere has a mass that linearly increases with its radius 
r and tends to oo as r cio. Clearly, such a system cannot represent a realistic DM halo of finite physical size. 
(2) The rotational speed in the isothermal sphere model increases linearly for small r, has an oscillatory dependence 
on r for intermediate r, and tends to a constant only asymptotically as r — >■ oo. Thus there is no strict reason 
why Vc,oo should be equal to Vc,q- (3) More importantly, the presence of the visible matter further complicates the 
situation; it contributes significantly to balancing the centrifugal forces of rotation at the solar circle. The density 
profile of an isothermal sphere is obtained as the solution to a second-order non-linear differential equation with 
appropriate boundary conditions at the origin. Thus the presence of a significant amount of visible matter within the 
solar circle will substantially change the nature of the solutions. (4) It is a common practice, in the context of the 
phenomenology of DM direct detection experiments, to truncate the speed distribution ([1]) at some chosen value of 
the local (solar neighborhood) escape speed Vescp', see, e.g., 0, ll])- However, this is not in general a self-consistent 
procedure because the resulting "truncated Maxwellian" speed distribution will not necessarily be a solution of the 
steady-state collisionless Boltzmann equation appropriate for a finite physical system. In addition, since the rotation 

curve for such a truncated Maxwellian will not in general be asymptotically flat, the relation (w^)^/^ = ^J^Vc^oo used 

to determine the value of (v^)^^^ in the Maxwellian speed distribution of the isothermal sphere, as done in the SHM, 
will not be valid. However, there are indeed well understood procedures that allow truncation of the DM halo at a 
finite radius and the velocity distribution at appropriate limits in a self-consistent manner. 

In this paper we construct a self-consistent model of the phase-space structure of the finite-size DM halo of the Galaxy, 
and study its implications for the analysis of the results of the direct detection experiments. The model is based on 
describing the phase-space DF of the DM particles in the halo by the so-called "lowered" (or truncated) isothermal 
models (often called "King models") [1^, which are proper self-consistent solutions of the collisionless Boltzmann 
equation representing nearly isothermal systems of finite physical size and mass. At every location x within the 
system a DM particle can have speeds up to a maximum speed i'max(x), which is self-consistently determined by the 
model itself. A particle of velocity Wmax(x) at x within the system can just reach its outer boundary, generally called 
the truncation radius, where the DM density by construction vanishes. Note that this maximum speed Wmax is not 
the same as the escape speed Ucscp, which is defined as the speed required for a particle to escape from the given 
location out to infinity. The escape speed is always larger than the maximum speed. For a King model DM halo of 
finite size, it is the maximum speed and not the escape speed that is relevant in considerations of direct detection 
experiments. As we shall see below, the speed distribution of the particles constitutin g a lowered isothermal model 
can be described as "quasi-Maxwellian" not unlike what is seen in recent simulations [25|. Apart from taking into 
account the finite size of the DM halo, we also explicitly account for the mutual gravitational interaction between the 
visible matter and the DM, as this interaction influences both the density profile and the velocity distribution of the 
dark matter particles. When the visible matter density is set to zero and the truncation radius is set to infinity, our 
halo model becomes identical to that of a single-component isothermal sphere used in the SHM. 



^ A recent study suggests a somewhat higher value of Vc,q ~ 250 km s ^ [26ll . which would imply a correspondingly higher value of 
(^"^y^^iBo ~ 306kms-i. 
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The model we develop here is along the lines described in earlier formulations [27|, |28| which included the effects of 
the gravitation of the visible matter on the dark matter in a self-consistent manner. Stated simply, these models 
account for the fact that the DM particles move in a gravitational potential to which both the visible matter and 
the dark matter make contributions. As the distribution of the visible matter may be fixed from observations, the 
parameters of the self-consistent velocity distribution function of the DM particles in presence of the visible matter 
can be determined by comparing the theoretically derived rotation curve with the observed rotation curve of the 
Galaxy. The velocity distribution function so determined is then used in the calculation of the direct detection rates. 

Of particular interest in this context is the rotation curve recently reported by Xue et al [2^ , which extends up to a 
galactocentric distance of ^ 60kpc. This rotation curve is derived from the kinematics of a sample of ~ 2400 blue 
horizontal-branch (BHB) stars taken from the Sloan Digital Sky Survey (SDSS) DR6 [sO] database. A noticeable 
feature of this rotation curve is that it gently falls from the adopted value at the Sun's location, v^.q ~ 220kms~^, 
to ^ 180kms~^ at ^ 60kpc. It is of considerable interest to study the implications of such a rotation curve for the 
direct detection experiments. In this paper, therefore, we use this new rotation curve for determining the parameters 
of our halo model. 

For the isothermal sphere model of the DM halo, the main effect of the gravitational infiuence of the visible matter 
on the structure of the DM halo is to shrink its 'core radius' [l^l, = (3(w^)/(47rGpo)) ^^^j where po is the central 
density ^ . This is because the DM is "pulled in" by the visible matter, increasing its effective central density. Thus 
the characteristic core radius of an isothermal DM distribution in presence of the visible matter is less than that of an 
isolated isothermal sphere whose properties are determined by its own gravity alone. The same, as we show below, is 
true for the lowered isothermal models. A related effect is the relative enhancement of the DM density near the plane 
of the disk due to additional gravitational influence of the disk. These effects are further discussed in the following 
section. 

We find that, for our best-fit halo model with parameter values /Odm,© = 0.2GcV/cm'^ and truncation radius rt = 
120 kpc, that provides the best fit to the rotation curve data of Ref. [1^ , the null result of the CDMS-II experiment , 
for example, gives a 90% C.L. upper limit on the WIMP-nucleon spin-independent (SI) cross section of ^ 5.3 x 10~® pb 
for a WIMP mass of ~ 71 GeV. 

We also study the issue of compatability of the claimed positive signal reported by the DAMA collaboration with 
the null results from the other experiments using, for illustration, the DAMA collaboration's original 2-bin annual 
modulation amplitude data. We find that for our best- fit halo model, there exists a region of the WIMP mass (m^) 
vs. WIMP-nucleon spin-independent (SI) cross section (cr^p = cr^Jj) parameter space bounded by 2.6 <, <, 10.8 GeV 
and 1.0 X 10~^ >, cr^p >, 1.2 x 10^^ pb, within which the DAMA's claimed modulation signal is compatible with the 
null results of other experiments. Similar "DAMA-compatible" regions of small WIMP masses obtain also for spin- 
dependent interactions. We, however, recognize that a more rigorous analysis using DAMA bins over smaller intervals 
should be able to better constrain the DAMA region within the context of our halo model, a task that we will take 
up in a subsequent work. 

The rest of the paper is arranged as follows: In section |ll] we describe our model of the dark halo of the Galaxy that 
includes in a self-consistent manner the gravitational effect of the observed visible matter on the structure of dark 
matter halo. The parameters of the model are determined by fitting the theoretically calculated rotation curves to a 
recently determined extended rotation curve of the Galaxy up to 60 kpc. The shapes of the dark matter density 
and mass profiles and the velocity distribution function of the DM particles are also derived. The implications of our 
model for the analysis of the results of WIMP direct detection experiments including the question of compatibility 
of DAMA's claimed positive signal with the null results of other experiments, are studied in section Hill Section ITVl 



^ The gravitational effect of the Galactic disk, of course, flattens the isothermal 'sphere'. Our formalism in this paper automatically and 
self-consistently incorporates this effect. Below, unless otherwise specified, by core radius we refer to that of the density distribution in 
the plane of the Galactic disk. Note, however, that while the resulting density distribution of the halo becomes anisotropic, one can still 
assume, as we do in this paper, that the velocity distribution of the DM particles remains isotropic everywhere, although anisotropic 
(velocity) models are also possible l20t . The effects of anisotropic and also triaxial halo models on the event rates in direct detection 
experiments and the resulting nature of the exclusion limits on WIMP parameters have been studied earlier; see, e.g., Ref. [sill and 
references therein. 
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summarizes the main results of this paper and concludes. 

In this paper we restrict our attention to only elastic scattering of WIMPs from nuclei due to spin-independent (SI) 
as well as spin-dependent (SD) WIMP-nucleus interactions. 

II. THE SELF-CONSISTENT TRUNCATED ISOTHERMAL MODEL OF THE DARK MATTER HALO 

OF THE GALAXY 

A. The Basic Formalism 

The Maxwellian velocity distribution ([1]) used in the SHM is simply related to the phase-space DF of the isothermal 
sphere 0, 



/is(x,v) = fiE) = ^° 3 exp [-E/a^] , 

(27rcr2)2 



(2) 



with {v ) = 3(T and E = <I'(x) + , the total energy (per unit mass) of the DM particle, $(x) being the total 
gravitational potential with the boundary condition $(0) = 0, so that the density at any point x is 



pis(x) = J /is(x,v)d3v = pocxp [-$(x)/tT2] ^ 



(3) 



where po is the central density. The isothermal sphere described by DF ([2]) is infinite in extent with a divergent total 
mass, and is characterized by two parameters, a and po- 

To describe "nearly isothermal" systems of finite size and finite total mass, one must have, in addition to the above two 
parameters of the isothermal sphere, a parameter that characterizes the finite size of the system. This is accomplished 
in King models by taking the DF to be of the following form [lO] : 



/k(x,v) = /(£) = 



pi(27rcr2)-3/2 (^e^Z-y^- _ 1^ for £• > , 
for f < , 



where 



f(x)^C-(i^,2 + $(x)Uvi.(x)-i«2^ 



(4) 



(5) 



is the so-called "relative energy" and ^'(x) (= — $(x) +C) the "relative potential" [2^. The three (constant) param- 
eters of the system are pi (with dimension of density) , a (with dimension of velocity) and the new parameter C (with 
dimension of potential or squared velocity) which is related to the finite size of the system (see below). 



The density at any position x is obtained by integrating /k(x,v) over all velocities giving 



Pk(x) = 



(2^(t2) 



2n3/2 



Pi 



which satisfies the Poisson equation 



exp 



*(x) 



erf 



exp 



4*(x) 



1 



V2$(x) = -V2«'(x) = 47rGpK(x) . 



- 1 
2^'(x) 



(6) 
(7) 

(8) 



At any location x the maximum speed a particle of the system can have is 

?^max(x) = V2*(x) , 



(9) 
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at which the relative energy £ and, as a consequence, the DF vanish. 



The King models approach the isothermal sphere solution at small radii r — |x|, but have density profiles that fall 
off faster than that of the isothermal sphere at large radii (which makes the total mass finite). Indeed, the density 
vanishes at r = r^, the truncation radius where £" = 0, representing the outer edge of the system. This is ensured by 
choosing 

C = $(rt), so that «'(rt) -$(rt) +C = 0. (10) 

As evident from Equation (fTO|). the parameter C fixes the finite size (r^) of the system. Actually, the King models 
form a sequence, each model being parametrized by the value of the "concentration parameter" Cja^. It can be 
shown [2^ that the sequence of King models goes over into the isothermal sphere in the limit C ja^ — s- 00 . 



Note that the usual escape speed, Woscp, defined as Wescp(^) — 2 [<&(oo) — $(x)] is related to the maximum speed 
i'max(x) defined by equation ([9|) through the relation 

<,,(x) = 2vl'(x) = 2 [$(rO - $(x)] 

= z;Lp(x)-2GM(r,)/rt, (11) 

where M{rt) is the mass contained within rt, with"* QM(rt) jrt = [$(00) — $(rt)]. Note also that, unlike in the case 
of the (infinite) isothermal sphere, the parameter a in the King models is not equal to the one-dimensional velocity 
dispersion of the particles constituting the system; the latter can be calculated for the DF given above and is a function 
of r, vanishing at r ~ rt- In fact, unlike in the case of the Maxwellian DF for which the velocity dispersion is linearly 
related to the most probable speed of the particles of the system, there is no simple relation between the velocity 
dispersion and the most probable speed of the particles of the system described by a skewed velocity distribution such 
as that described by the King DF (|3]), and as such, the velocity dispersion has no special significance and does not 
uniquely specify the full DF of the King model. 



In the subsequent discussions we shall take the King model DF ^ as the DF of the DM particles constituting the 
finite-size DM halo of the Galaxy, in the rest frame of the Galaxy (the subscript K will be dropped for convenience) . 



For a "pure" DM halo of finite size represented by the DF ([4]) without including the gravitational influence of 
the visible matter embedded within the halo, the gravitational potential $ in the above equations can be written as 



<I>(x) = ^^^(r), where ^oM{r) is the spherically symmetric self-gravitational potential of the pure DM halo represented 
by the DF ([4]). In reality, however, a test particle sees the total gravitational potential of both the DM and the visible 
matter embedded within the DM halo. Thus, we should write 



$(x) = $oM(x) + $.i„(x), 



(12) 



where now $dm represents the contribution of the DM component to the total gravitational potential in presence of 
the visible matter, and is the gravitational potential of the observed visible matter of the Galaxy. These potentials 
satisfy their respective Poisson equations 



V2$om(x) = 47rGpoM(x) , and V^'^^.^'x) = AnGp^^^) , 



(13) 



where the DM density is 



Pdm(x) = pi 



exp 



4*(x) 



21'(x) 



(14) 



with 



*(x) = [$om(|x| =rt)-^$_(|x| 



$Om(x) +$vis(x) 



(15) 



Note that wc have chosen the boundary condition "J?{0) = 0, so that <I>(oo) is a non-zero positive constant. 
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using (|10p . and p^^^ is the known visible matter density of the Galaxy. 



Note that in presence of the visible matter, whose spatial distribution is non-spherically symmetric, the DM spatial 
distribution and hence its gravitational potential will also be non-spherically symmetric in general. However, for large 
enough values of the truncation radius (jt) of the halo (more specifically, for large compared to the visible matter 
disk's scale length — the latter being ^ few kpc; see below), both $om and <i>via become close to spherically symmetric 
at |x| — > rt, which we also see in our numerical calculations described below. On the other hand, as already mentioned 
in footnote 3 in Section U the velocity distribution will be assumed to remain isotropic. 



The visible matter distribution p^^^lx) (and hence the potential $^ig(x)) being known from various observational data 
and modeling, solutions of equations (|13p . (jl4p . (|15p with appropriate boundary conditions, which we choose as 

$dm(0) = $™(0) = 0, and (V$oM)|x|=o = (V*vis)|,|^o = 0' (16) 
give us a three-parameter family of self-consistent solutions for jOom(x) and $om(x) for chosen values of the parameters 



For a given solution, the rotation curve of the Galaxy, Vc{R), i.e., the circular rotation velocities in the equatorial 
plane of the Galaxy as a function of the Galactoccntric distance (i?) on the equatorial plane, can be calculated from 



$total(i?,2 = 0) 



$oM(i?,2 = 0) + $_(ii',z = 0) 



(17) 



z being the distance normal to the equatorial plane. The best-fit values of the relevant parameters can then be 
determined by comparing the theoretically calculated rotation curves with the observed data on the rotation curve of 
the Galaxy. The parameter values so determined are then to be used as inputs to the analysis of the results of the 
direct detection experiments. 



We have developed an efficient numerical scheme to solve the coupled non-linear Poisson equation (fT3|) for DM with 
Pom given by equations ()14p and ()15|) with boundary conditions ()16|) , for a given potential-density pair of the visible 
matter, through an iterative procedure first discussed in Each stage of this iteration procedure involves a Poisson 
solver which is tested on a number of exact analytical formulae for potential-density pairs given, for example, in [20| . 
The iteration process converges typically within ten iterations. We have also verified that this iterative scheme yields 
good agreement with the known results for both the non-singular isothermal sphere and King Models given in [20| . 
to within a few percent. 



We assume that the density distribution of the visible matter can be effectively described by the well-known model of 
a spheroidal bulge superposed on an axisymmetric disk [s^ . The density distributions of these components are given, 
respectively, by 

/ 2 \ -3/2 

p,{r)^p,{0)i^l + ^j , (18) 

and 

p,(,) = ||e-(«-«")/«'' e-l^l/\ (19) 

wherer = (i?2 + z2)i/2 Typical parameter values are HHH] ps(0) = 4.2 x 10^ Mg pc"^ , a = 0.103 kpc , i?d = 3 kpc , 
and h = 0.3 kpc, with Rq — 8.5 kpc, the solar Galactocentric distance, and Sq 48M0pc~^, the surface density 
of the disk at the solar location. The expressions for the gravitational potentials, (j)^ and 0^, corresponding to above 
forms of ps and pd, are given in [3^ . Alternatively, they can be directly calculated using our numerical Poisson solver. 
This specifies the ^^^^ — (ps + (pd used in the numerical calculations in this paper. 



There exist more detailed models of the mass distribution of the visible matter in the Galaxy than those adopted 
above, involving a stellar "thick disk" (with a scale height of ~ 1 kpc) and the disk formed by the Galaxy's interstellar 
medium (ISM) in addition to the stellar spheroid and the "thin" disk described by equations ([T5)) and (|19p . However, 
there are also large uncertainties in the values of the parameters that characterize these various components, as clear 
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from the values of the parameters of the two representative models summarized, for example, in Table 2.3 of Ref. [20| . 
We find, as shown below, that the typical values of the various visible matter parameters chosen above give reasonably 
good fit to the rotation curve data in the inner Galaxy region, R <, Rq w 8.5 kpc, where the effect of the Dark Matter 
should be minimal. As such, we believe our "minimal" model of the distribution of the normal matter of the Galaxy 
adopted above is good enough for the purpose of illustrating the general nature of effects arising from a self-consistent 
description of the Galaxy's phase space structure that includes the mutual gravitational interaction of the visible 
matter and the Dark Matter components of the Galaxy. The formalism of this paper and the numerical procedure 
adopted by us are, however, quite general and applicable to any given model of the mass distribution of the Galaxy. 

In our numerical calculations, we have taken the DM density at the solar location, Pdm.g = Pr>m{R = -RojO), as the 
"observable" density parameter of the King model instead of the parameter pi in equation ([T4| . We have run a 
large number of numerical models of the Galaxy for three different values of Pdm,© = 0.2 ,0.3 and 0.4 GeV/ cm'^ , in 
each case with a wide range of values of the other two King model parameters rt and a and with the visible matter 
parameters fixed as described above, to determine the parameters that best describe the Galaxy's rotation curve data 
mentioned above. 



B. Dark Matter Density Profile 

The main effect of the gravitational influence of the visible matter on the DM halo, namely, increased central concen- 
tration and reduced core radius of the DM density proflle, is illustrated in Figured] The left panel of Figure [T] shows 




R (kpc) R (kpc) 

FIG. 1: The density profiles of the lowered (truncated) isothermal dark matter halo described by the King model DF, eq. ([4|, 
(a) in absence ("uncoupled"), and (b) in presence ("coupled"), of the visible matter described by equations (|18|) and (|19[) . 
for three difi'erent sets of values of the King model parameters (rj, cr, Pdm,©) that yield good fits (in the coupled case) to the 
rotation curve data of the Galaxy up to ~ 60 kpc. In (b) the density profiles refer to those on the equatorial plane (z = 0) of the 
Galaxy. The inset in (b) shows the ratio of the density on the z = Q plane to that on the z axis as a function of galactocentric 
distance in the coupled case, for the same chosen sets of the King model parameters as indicated. 

the density profiles of a 'pure' spherically symmetric finite-size DM halo described by the King model DF ^ without 
including the gravitational influence of the visible matter on the DM (hereafter referred to as the "uncoupled" case), 
whereas the right panel shows the DM density profiles (in the plane of the Galactic disk (z = 0)) under the additional 
gravitational influence of the visible matter (hereafter referred to as the "coupled" case). The curves in both panels 
are for three different sets of values of the King model parameters {rt, a, Pum.q) as indicated, that yield good fits (in 
the coupled case) to the rotation curve data of the Galaxy up to ^ 60 kpc (see discussions below and Figure [2|) . The 
increased central density and the reduced core size of the density profiles in the coupled case relative to those in the 
uncoupled case are clearly seen. 

Another direct effect of the gravitational influence of the visible matter on the DM is the enhancement of the DM 
density on the plane of the Galactic disk (z = 0) relative to that off the disk. The inset in the right panel of Figure [T] 
shows this effect clearly where it is seen that this disk enhancement of DM density can be as large as ~ 30 - 40 %. 
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The typical scale length of this "dark matter disk" is a few kpc. 

C. Rotation Curves 

Figure [2] shows our theoretically calculated rotation curves for the Galaxy with its DM halo described by the King 
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FIG. 2: Rotation curves for the Galaxy with its DM halo described by the King model DF, eq. (|4]) and including the gravitational 
effect of the visible matter described by equations (|18p and (|19p . for three different sets of values of the King model parameters 
(n , (7 , Pdm.o) as indicated. The curves are shown for three different values of Pdm,©, and in each case, the curve shown 
corresponds to the values of the other two parameters (rt and cr) that yield best fit to the rotation curve data of the Galaxy 
estimated in Ref. [29| . The data with error bars are from Table 3 of Ref. [2^ (their data set Kir.ii corresponding to their 
simulation II) for galactocentric distances R > 7.5 kpc, and from Ref. (3^ for R < 7.5 kpc. 



model DF, eq. ^ and including the gravitational effect of the visible matter described by equations and ([T^. for 
three different sets of values of the King model parameters (rt, cr, poM.o) as indicated. For comparison, the rotation 
curve data for R up to 55 kpc given, for R > 7.5 kpc, in Table 3 of Ref. [2^ (their data set Vdrji corresponding to 
their simulation II) arc shown together with the data from Ref. [U for R < 7.5 kpc. The theoretical rotation curves 
are shown for Pdm.o — 0.2 , 0.3 and 0.4 GeV/ cm'^, and in each case, the shown theoretical curve corresponds to the 
values of the other two parameters (rt and a) that yield best fit (giving lowest x^) to the above mentioned rotation 
curve data of the Galaxy. 



For the range of King model density parameter considered (0.2 < Pdm,© < 0.4 GeV/ cm'^), the "global" best fit to 
the rotation curve data (giving globally lowest value of x^) is obtained for the King model DM parameter values 
Pdm.0 = 0.2 GeV/ cm^, rt ~ 120 kpc and a ~ 300kms~^ (the dashed curve in Figurc[5|) ^. The separate contributions 
of the visible matter (VM) and the Dark Matter (DM) components to the total rotation curve and the total mass of 
the Galaxy for this "best- fit" set of parameter values are shown in Figure [3] 



It is of interest to note that the DM halo mass for the best-fit model obtained above. A/dm ^ 1-3 x 10^"'^ M0, is 
only marginally more than the total visible matter mass, AfvM ^ 8.4 x 10^° A/©, giving a total Galaxy mass of 



^ Note that, as already mentioned, the parameter cr in the King model is not to be confused with the one-dimcnsional velocity dispersion 
of the DM particles. 
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FIG. 3: The contributions of the visible matter (VM) and the dark matter (DM) components to (a) the total rotation curve (left 
panel) and (b) total mass (right panel) of the Galaxy for the best-fit King model DM halo parameters Pdm,© = 0.2GeV/cm'^, 
rt = 120 kpc and a ~ 300kms~^. The rotation curve data in the left panel are the same as in Figure (2] 



~ 2.1 X 10^^ Mq for the adopted VM model. Actually, this is fairly independent of the exact values of the VM model 
parameters adopted as long as those parameters are such that the VM by itself gives the dominant contribution to the 
rotation curve data at small R Rq . The main reason for the relatively low DM halo mass is the declining nature 
of the rotation curve beyond the solar circle. Within the context of the truncated isothermal models of the DM halo 
described by the King models studied here, higher DM halo mass models — which can be obtained by choosing larger 
values of one or more of the parameters Pdm.oi ft or cr than those of the best-fit model obtained above — typically 
lead to rising or non-declining rotation curves beyond the solar circle, and are thus unable to explain the declining 
rotation curve data. 



D. Maximum Speed and Speed Distribution of the particles 

The self-consistcntly calculated maximum speeds of DM particles in the halo as a function of R arc shown in Figure [4] 
(left panel) for the same three sets of King model halo parameters as in Figure [2l The corresponding normalized speed 
distribution functions of the DM particles (in the Galactic rest frame), f{v) = ^^/(x, v) (with J f{v)dv = 1), at 

Sun's location (i? = i?0, z = 0), are also shown (right panel). For comparison, the Maxwellian speed distribution for 
the Standard Halo Model (SHM) at sun's location is also shown. Evidently, the speed distribution for the truncated 
"isothermal" halo is significantly non-Maxwellian, especially at the high speed end of the distribution, not unlike the 
behavior found in recent numerical simulations |25| . 

For the three King model parameter sets considered in the left panel of Figure |4] the maximum speeds of the DM 
particles at the location of Sun are, from bottom to top, 396, 403 and 439 km/s, respectively. These values and the 
speed distribution functions shown in Figure |4] will be used below in the calculation of the expected WIMP detection 
rates and analysis of the results of the direct detection experiments. We emphasize that the above values refer to the 
maximum speeds that a dark matter particle can have at the location of the Sun. Since stars in the Galaxy do not in 
general follow the same velocity distribution function as the dark matter particles, they can have speeds larger than 
the maximum speeds of the dark matter particles mentioned above. Stars with speeds larger than the escape speeds 
determined by the total mass and radius of the halo (see equation (fTTj) ). will, of course, eventually escape from the 
Galaxy. 
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FIG. 4: Maximum speed Wmax(x) = ^/2W(x) in the Galactic equatorial plane as a function of Galactocentric distance R (left 
panel) and the normalized speed {v — |v|) distribution function of the DM particles, f{v) = ^^^/(x, v) (right panel), at Sun's 
location {R = Rq, z = 0), for the same three sets of King model parameters as in Figure [5] as well as for the Maxwellian speed 
distribution in the case of Standard Halo Model (SHM). 



III. IMPLICATIONS OF THE SELF-CONSISTENT TRUNCATED ISOTHERMAL HALO MODEL FOR 
THE ANALYSIS OF WIMP DIRECT DETECTION EXPERIMENTS 



A. Direct Detection of WIMPs: The Basic Formalism 



For simplicity we restrict our attention in this paper to the case of elastic scattering of the WIMPs from nuclei. 
Consider a WIMP (hereafter denoted by x) of mass scattering elastically from a nucleus of mass M. The nucleus 
receives a recoil energy, = q^/{2M) = {ij?v? /M){1 — cosO*), where 6* is the scattering angle in the centre of 
momentum frame, q the momentum transferred to the recoiling nucleus, /i = m^M / {m^ + M) is the WIMP-nucleus 
reduced mass and u = |u| is the speed of the WIMP relative to the nucleus assumed to be at rest on earth. The 
minimum WIMP speed that can produce a recoil energy of the nucleus is 

UnuniEn) = (^^) ^ • (20) 
The differential event rate per unit detector mass (typically measured in counts/day/kg/keV) can be written as @, |1] 

-^(i?^,t) = -^p,e(£;«,t), (21) 

where = Pdm,o is the density of WIMPs in the solar neighborhood, (T{q) is the WIMP-nucleus effective interaction 
cross section, and 

aER.t)^ [ d'ul^. (22) 

v'u>u,„i„(£;fj) 

Here /(u, t) is the time-dependent velocity distribution of the WIMPs in the solar neighborhood relative to the detector 
at rest on earth, and Umin is defined in equation (PU)) . Recall that /k (x, v) defined in equation (U) refers to the rest 
frame of the Galaxy, and has no explicit dependence on time. The time dependence of / arises from the motion of 
the earth with respect to the Galactic rest frame. The two distribution functions are simply related as 

/(u,i) = — /K(x = xo,v = u + VE(i)) , (23) 
Px 

where Xq represents the sun's position {R = 8.5kpc, z = 0) and v^it) is the earth's velocity vector in the Galaxy's 
rest frame. Note that the DF /k vanishes for speed v > Umax defined in equation ^ . 
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Using equation (123)) we can explicitly write equation (|22|) as 

<.U„ax(cOSe) 

^{Eii^t) / ^(cOSe) e (Uinax - Umin) / u/k (x = Xq , V = U + Ve (t ) ) , (24) 

where Mmax is the positive root of the equation 

■(^max = «iax + «E + Sw^axWE COS 6* , (25) 

with Vmax given by equation ([9]) evaluated at x = X0, and ue = |ve| given by [l^ 

VE{t) W0 + i;E,orbC0S7C0S [w(t - to)] , (26) 

where Vq « fc,0 + 12kms~^ is the magnitude of the dominant component of sun's total velocity vector (which includes 
its circular velocity around the Galaxy of magnitude Wc,o ~ 220 km s^^ and its "peculiar" velocity of magnitude 
« 12kms~^) in the Galactic rest frame, we, orb ~ 30kms~^ is the average orbital speed of the earth around the sun, 
cos 7 ss 0.51, w = 27r/365 and to = 152.5 (corresponding to 2nd June when the earth's speed with respect to the 
Galactic rest frame is maximum), t being counted in days. It is the periodic variation of ve, equation (|26|1 . that gives 
rise to the annual modulation signal claimed to have been detected by the DAMA collaboration. 



For detectors composed of more than one kind of elements, the total differential event (i.e., nuclear recoil) rate is 



(^) -TfJ4^) , (27) 



where fi is the mass fraction of, and (^j^^ _ the differential event rate (eq. (PT|) ) for, the nuclear species i, with 
corresponding nuclear mass Mi, WIMP-nucleus cross section <Ji{q) and WIMP-nuclcus reduced mass fii. 



The number of nuclear recoil events in a recoil energy range between Ej^ and is 



NniElEl) = Yl fj dEn (^^) IC.{En) 



(28) 



where ICi — AiiTe{Eji) is the total exposure of the detector, Aii being the total mass of the species i in the detector, 
T the total exposure time, and e{Eji) the energy-dependent detector efficiency. 

The WIMP-nucleus total effective scattering cross section, a{q), can be written as a sum of two contributions arising 
from spin- independent (SI) and spin-dependent (SD) effective couplings of the WIMP to the detector nucleus [2l,[l,[l^ : 

friq) ^ asi{q) + asu{q) ■ (29) 



Spin-independent (SI) (or Coherent) Scattering : In this case, it is assumed that the WIMP interacts coherently with 
the nucleus as a whole. The cross section is generally written as 0| 

aiq) ^ ao\F{q)\^ (30) 

where cto is the zero-momentum WIMP-nucleus scattering cross section, and F(q) is a momentum dependent form 
factor that arises from the finite size of the nucleus. The form factor F(q) is normalized to i^(0) = 1. The coherent 
nature of the interaction is implemented under the usual assumption that ctq scales with the square of the atomic 
mass number A (number of protons plus neutrons) of the nucleus. For purely scalar interactions, with equal WIMP 
couplings to protons and neutrons, one can write (Jq in terms of the WIMP-proton (or WIMP-neutron) effective cross 
section a^p = a^n as 

-o'=4\(-^X^'^ (31) 
where fi^p is the WIMP - proton reduced mass. 



13 



Wc shall use the conventional Helm form of the nuclear form factor 0, [l^ : 

N o -q^s'nt? sm{qro/h) - (grp/h) cos{qro/h) 

with s = 0.9 fm. the nuclear skin thickness and tq ~ 1.14^^/'^ fm, the effective nuclear radius. 



(32) 



Spin- dependent (SD) Scattering : In this case the WIMP couples to the total spin J of the nucleus which has 
contributions from the spins of the individual protons and neutrons within the nucleus. Following Ref. [Toj we use 
the effective WIMP-nucleus SD cross section generically written in the form [35 1 



o'SD(g) 



32^1, 
2J + 



(33) 



where Up and a„ are respectively the axial four-fermion WIMP-proton and WIMP-neutron couplings in units of 
2-\/2G_F [36l-l38l . The nuclear structure functions Spp{q), Snn{q) and S^ n (q) are taken, for Ge from Ref. [3§|, for Si 
from Ref. [401 j ^r Al from Ref. 41| , and for Na, I and Xe from Ref. [i^ . 



In this paper, following the standard practice, we shall study the situations when either On = or = ; the relevant 
effective cross section in the two cases will be denoted by and cr^, respectively. In general both a,i and Op can 
be non-zero, but we will not consider this case here. 



B. Results of Direct Detection Experiments and Analysis Techniques 



To illustrate the general implications of our halo model for the WIMP mass and cross section as implied by the results 
of direct-detection experiments, wc consider in this paper the results of three of the "null" experiments (i.e., those 
without any definite claim of detection so far)_^amely, CDMS-II [ill,!!!, CRESST [13 and XENONIO [ll,[ll, in 
addition to those of the DAMA collaboration [10| which has claimed a positive signal based on the claimed detection 
of the annual modulation of the recoil event rate. As we shall see, results of these three null experiments together 
provide useful information regarding the question of compatibility of the positive results of DAMA with the negative 
results of other experiments over a wide range of the WIMP mass from ^ 1 GeV to 100 GeV. Below, we first 
describe the DAMA results and the analysis procedure wc follow in order to derive the constraints on WIMP mass 
and interaction cross section imposed by the results, and then discuss the same for the null experiments. 



DAMA/Nal and DAM A/LIBRA : 

The DAMA coUaboration has claimed detecting a non-zero annual modulation signal at a confidence level (C.L) 

of 8.2(7 in their number of detected nuclear recoil events over a period of seven (for the DAMA/Nal experiment [l]), 
plus four (for the DAMA/LIBRA experiment [l3|), annual cycles, with a total exposure of 0.82 ton-year^. The 
DAMA collaboration attributes this annual modulation to the periodic variation of we, equation ([26]), which results 
in a periodic variation of the recoil event rate (see discussions in the previous subsection). 



Given the time dependence of we , equation (|26p , the time dependence of the differential event rate (|21[) due to motion 
of the earth around the sun can be approximately written as [^, [3, Il9l | 



dn 

dE^. 



{ER,t) « So{Er) + Sm{ER) cos oj{t^ to), 



(34) 



where Sq is the average recoil rate over a year and Sm is the "modulation amplitude" defined as 



Sm{ER) - - 



'"^ {Er, June2) - {Er, Dcc2) 



dER 



dER 



(35) 



While this paper was being prepared, the DAMA/LIBRA collaboration released their latest results with an additional exposure of 0.34 
ton-year corresponding to two additional annual cycles [44| . With this, the statistical significance of the annual modulation signal is 
now 8.9cr for the cumulative exposure. Wo do not include these new results in this paper. 
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A nonzero value of Sm is taken to be a signal for WIMP-induced recoils. 

To compare with DAMA data given in specific recoil energy bins, the average value of Sm over a given energy range 
is calculated as: 

^2 — i^l J El 

A scintillation detector such as that used in the DAMA experiment directly detects only the part of the recoil nucleus 
energy that goes into electromagnetic channel and produces the observed scintillation. The actual energy of the recoil 
nucleus Eji is related to the detectable energy from the scintillation light yield, E^, — often called the "electron- 
equivalent (ee) energy" and denoted by "keVee" with energy in units of keV — through the relation Ed = QEr, 
where the "quenching factor" Q{< 1) depends on the nuclear material composing the scintillation detector. For the 
DAMA detector (composed of Nal crystal), ~ 0.3 and Q\ « 0.09. At the same time, as first pointed out in 
Ref. [i^ and studied in detail in the context of the DAMA detector in Ref . [1^ , for certain energies and incidence 
angles of the particle (for example, along the crystal axis) , the incident particle transfers energy only to the electrons 
(rather than to the nuclei) of the scintillator material. In such a case, called "channeling", one has Q ~ \. We shall 
use the following simple parametrizations, suggested in Ref. [13], of the fractions of "channeled" events for recoiling 
Sodium and Iodine nuclei, based on simulation results given in Ref. [46| for the DAMA experiment: 

^ l + 1.14Sfl(kcV) ' ^ l + 0.75^fl(keV) ' ^^^^ 

Note that the use of keVee is only a bookkeeping device to distinguish the actually measured energy by the detector 
from the true recoil energy, the actual energy unit still being keV. For detectors for which the concept of channeling 
does not exist (for example, CDMS or XENON), the quenching factor can be included in the energy calibration and 
final data quoted directly in terms of the recoil energy. 

For the analysis of the DAMA results, in this paper we consider for simplicity only the 2-bin data set given by the 
DAMA collaboration [l3|, namely, the low-energy bin 2 ~ 6 keVee within which a non-zero modulation amplitude is 
measured, and the high energy bin 6-14 keVee in which the modulation amplitude is consistent with zero. Table I 
gives the modulation amplitudes [Sm] measured by DAMA in these two energy bins. 



Energy (keVee) 


Modulation Amplitude {Sm) (counts/day/kg/keVee) 


2-6 
6-14 


0.0131 ± 0.0016 
0.0009 ± 0.0011 



TABLE I: DAMA modulation amplitude data [K 



The efficiency of the DAMA detector is taken to be unity. To compare with the actual experimental data we take into 
account the finite energy resolution of the detector whereby the actually measured energies are taken to be normally 
distributed about the true detectable energy Ejj (i.e., the energy that would be measured if the detector had 100% 
energy resolution) with a standard deviation (t{Ed) = (0.448 keV)-\/ED7keV -|- 0.0091£'£) . Thus, the expected 
modulation amplitude is obtained by convolving equation p5|) (after changing the integration variable to Ed ~ QEr) 
with a normalized Gaussian with the above standard deviation. The expected modulation amplitude over a given 
interval of measured energy between Edi and Ed2 is then calculated from equation p6p after the above convolution. 

With the theoretically expected modulation amplitude in the fc-th energy bin, 5*^^,, calculated as described above, 
we perform, following the simple analysis procedure of Ref. (isj , a fit to the experimental modulation amplitude 
data given in Table I with 

^^^y(^^ii±z^]\ (38) 
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where Sm,k is the experimentally measured modulation amplitude in the fc-th energy bin and ak the corresponding 
error in the experimental value given in Table I. In this simple analysis procedure, there are only two free parameters 
in the problem, namely, the WIMP mass and the relevant WIMP-nucleon cross section a^. The latter is cr^p in 
the SI case (see eq. (|3ip ). and either cr^^ or o'^^ in the SD case (see eq. (|33p and the discussions following it). For a 
given value of we find Xmim minimum value of x^, by scanning over the values of the relevant cross section 
a^. The 90% C.L. allowed region of the relevant cross section, for the given value of m-^, is then found by accepting 
those values of a-)^^ for which "LXnin — 2.71, provided Xmin ^ 2. More rigorous multi-parameter analysis procedures 
described, for instance, in Ref. |19j . may alter the constraints derived here at some quantitative level. However, we 
believe the simple analysis procedure adopted here is sufHcient for our main purpose of illustrating, in the context 
of our self-consistent model of the dark halo of the Galaxy, the general nature of the constraints in the (m^, tr^) 
parameter space implied by the experimental results. The resulting constraints imposed by the DAMA results are 
described in the next subsection together with those imposed by the null results of the other experiments. 

Experiments giving null results : 

Table II summarizes the relevant features of the "null" experiments we consider in this analysis. Details of each 
experiment can be found in the cited References. 



Experiment 


Target 


EfTective exposure (kg-days) 


Threshold (keV) 


CDMS-II [llj 


Ge 


304.5 (SI), 23.5 (SD) 


10 


CDMS-II [43| 


Si 


12.1 (SI), 0.57 (SD) 


7 


CRESST-I [12] 


AI2O3 


1.51 (SI), 0.53 (SD) 


0.6 


XENONIO [13, JJ 


Xe 


136 (SI), 64.7 (SD) 


6.1 



TABLE II: Relevant features of the experiments with null results considered in this paper. 



For the purpose of analysis, the null experiments can be divided into two classes: (a) Those which report no events at 
all (after all the relevant cuts are employed) , and (b) those which report events (after analysis cuts) but ascribe them 
to background. In both cases, we follow the simple analysis procedure outlined in Ref. For the former case (a), 
we derive 90% C.L. upper limit on the relevant WIMP cross section cr,^ for a given WIMP mass by using Poisson 
statistics and demanding that the theoretically predicted number of events, N, over the entire energy range of interest 
be such as to allow a Poisson probability of observing zero events as bad as 10% but not worse. This corresponds 
to Nim-y, Q-y) = 2.3. In the case (b), a simplified version of Yellin's [4§| optimum interval method, as suggested in 
ReL [l8[, is used, whereby we again calculate, for a given mass m^, and for all contiguous combination of energy bins, 
the 90% C.L. upper limit allowed theoretically expected number of events N{m^,a^) corresponding to the observed 
number of events n as dictated by Poisson statistics, i.e., by the formula X]j=o "T" ^^P(~-^) = ^.1, and choose the 
most stringent (i.e., the lowest) upper limit value of so obtained. 

CDMS-II (Ge) [m : 

This is currently the experiment with the largest effective exposure amongst all currently running experiments. Two 
events were observed, after applying all cuts, in the signal region at recoil energies 12.3 keV and 15.5 keV in the 
10 - 100 keV window. The effective quenching factor is already included in the energy calibration. Since the 
estimated probability of observing 2 or more background events in this window is ^ 23%, these events were ascribed 
to background. We obtain upper limits on the relevant WIMP cross section as a function of WIMP mass using the 
simplified version of Yellin's |49| optimum interval method described above. In this analysis we have used a minimum 
bin width of 0.4 keV. The total effective exposure used in this calculation includes the spectrum-averaged equivalent 
exposure of w 194.1 kg-days (for a WIMP mass of ^ 60 GeV/c^) reported in Ref. [ll| combined with that of CDMS-II 
collaboration's previous published paper [50|, ~ 110.4 kg-days (after reduction of the figure quoted in Ref. Q by a 
factor of ~ 9% due to improved estimate of their detector mass as prescribed in Ref. [ii|). The above figures include 
the approximately constant (in energy) detector efficiency of ^ 30% already folded in. 

For spin-dependent interaction, only the ^"^Gc isotope, whose natural abundance is w 7.73%, is sensitive to the 
WIMPs. The effective exposure for SD interaction is correspondingly reduced and is significantly less than that for 
the SI interaction case. 
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CDMS-II (Si) El : 

This experiment employs a high purity sihcon crystal as the detector material. No events passed all the data analysis 
cuts. We derive 90% C.L. upper limit on the relevant cross section as a function of the WIMP mass by using Poisson 
statistics discussed above, i.e., by demanding that, for a given WIMP mass jti^, the upper limit value of cross section 
yield 2.3 events. 

For SD interactions, only the ^^Si isotope, whose natural abundance is 4.68%, is effective. 
CRESST-I [13 : 

The Phase I of the CRESST experiment employed sapphire (AI2O3) detectors with an exposure of 1.51 kg-days over 
an energy range of 0.6-20 keV. Measured energies are calibrated to recoil energy ((5 = 1). All the observed events were 
ascribed to background. We take into account the energy resolution of the detector descri bed by a Gaussian distribu- 
tion of the measured energies around the true value with a standard deviation (j{Eb) = \/ (0.220 kcV)^ + (0.017i?ij)^) 
[T^^. U ppe r limits were obtained similar to CDMS-II(Ge) analysis with a smallest allowed energy-interval width of 
1.2 keV [ill. 

For SD case only the Aluminum is sensitive since Oxygen has J=0. 
XENON 10^^ : 

XENONIO is a dual phase (liquid and gas) xenon time projection chamber and uses the ratio of ionization to scintilla- 
tion yield for discriminating between the dominant electron-recoil background and the looked for nuclear-recoil WIMP 
signal over a nuclear-recoil energy range of ^ 6.1-36.5 keV (Q — 1) with the relevant energy calibration parameter 
^eff = 0.14 We use an energy resolution oi a{ER) = {0.579keY)^/E]^JkeV + 0.021Er as suggested in Ref. 
A total of 10 events (without background subtraction) were recorded in the WIMP signal region for a total effective 
exposure of about 136 kg-days after analysis cuts. As in the case of the CDMS-II Ge analysis described above, we 
again obtain upper limits on the relevant WIMP cross section as a function of WIMP mass using the simplified version 
of Yellin's optimum interval method described above. 

For spin-dependent WIMP interaction, only the isotopes ^^^Xe (spin-1/2) with a natural abundance of ^ 26.4% and 
"'^^"'^Xe (spin-3/2) with a natural abundance of ~ 21.2% are effective. 

C. Implications for the WIMP parameters : Exclusion plots and compatability of null results with DAMA 

signal 

Figures [S] - [7] show our main results in terms of the constraints on the relevant WIMP-nucleon cross sections as a 
function of the WIMP mass, as implied by the results of the direct detection experiments. In all these Figures, the 
regions above the curves are excluded at the 90% C.L. by the results of the respective experiments. 

Figure [5] shows our 90% C.L. upper limits on WIMP-nuclcon spin-independent cross section as a function of WIMP 
mass as implied by the latest CDMS-II results |ll| . For our halo model with Pdm.o = 0.3 GeV/ cm^ (with other halo 
parameters fixed by fitting to the rotation curve data) , the lowest upper limit on the WIMP-nucleon spin- independent 
cross section is 3.4 x 10~* pb at = 69 GeV — coincidentally almost identical to the value of 3.8 x 10~^ pb quoted 
in Ref. [ll| which uses SHM in their analysis (presumably with a chosen value of the Galactic escape speed of 
600kms~^ [1]). However, for our self-consistent halo model that gives the best fit to the new rotation curve data — 
this model has Pdm,© = 0.2 GeV/ cm^ — the above upper limit on the SI cross section changes to ^ 5.3 x 10~* pb at 
— 71 GeV. The upper limit value of the cross section approximately scales inversely with the value of Pdm,o- 

Figures m and [7] show the results of our study of the compatibility of the DAMA experiment's claimed positive result 
(based on their original 2-bin annual modulation amplitude data) with the null results of various other experiments, for 
the cases of WIMP-nucleon spin-independent and spin-dependent cross sections, respectively. For the spin-independent 



^ Ref. Il2)| actually quotes their energy resolution in terms of the FWHM, AEpi^'^^^ , which is related to the standard deviation a by 
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FIG. 5: 90% C.L. upper limits on the WIMP-nucleon spin-independent cross section as a function of WIMP mass as implied 
by the latest CDMS results [Til], foi' our self-consistent truncated isothermal (King) model of Milky Way's DM halo described 
by the same three sets of King model parameters as in Figure [2l The WIMP speed distribution function and the maximum 
WIMP speeds for these King model parameters used in deriving the above upper limits are as shown in Figure |4l 



case, we find that there is a range of small WIMP masses, 2.6 ^ uiy, ^ 10.8 GcV, within which DAMA's claimed 
modulation signal is consistent with the null results of other experiments. The allowed WIMP-nucleon SI cross 
section varies from ~ 1.0 x 10~^ pb at the lower end of the WIMP mass range to ^ 1.2 x 10^^ pb at the upper mass 
end. For the spin-dependent case, the DAMA-compatible mass range for WIMP-neutron interaction (i.e., Op = 0) is 
2.5 <, rriy^ <, 7.4 GeV, with the corresponding WIMP-neutron interaction cross section range of 9.5 x 10"^ - 42.0 pb. 
For the case of WIMP-proton interaction (with a„ = 0), the corresponding WIMP mass and cross section ranges are 
2.3 <m^< 16 GeV and 79.0 - 0.45 pb, respectively. 

We emphasize that the DAMA-compatible ranges of WIMP mass and cross sections obtained above are based on our 
simple analysis procedure using the original 2-bin DAMA annual modulation amplitude data, and are comparable to 
those obtained earlier within the context of the SHM, e.g., in Ref. [19| ( their 2-bin analysis results). More rigorous 
analyses using DAMA bins over smaller intervals, as done, e.g., in [1^, should be able to place more restrictive 
constraints on the DAMA-compatible ranges of WIMP mass and cross section within the context of our halo model. 
This will be taken up in a work now in progress. 



IV. SUMMARY AND CONCLUSIONS 



The standard halo model (SHM) of the dark matter halo of the Galaxy, customarily used in the analysis of the results 
of WIMP direct detection experiments, is not suitable for describing a finite-size system such as the Galaxy. In the 
SHM, the DM halo of the Galaxy is described as an isothermal sphere with a MaxwcUian velocity distribution of its 
constituent particles. The isothermal sphere solution of the steady-state coUisionless Boltzmann equation is actually 
infinite in extent and has a divergent total mass with the mass contained inside a radius r increasing linearly with 
r [2^. In the literature on direct detection of WIMP dark matter particles, it is a common practice to truncate the 
Maxwellian speed distribution assumed in the SHM at a chosen value of the local escape speed of the Galaxy. This, 
however, is not a self-consistent procedure because the resulting truncated speed distribution does not in general 
satisfy the steady-state coUisionless Boltzmann equation which a system of coUisionless particles such as WIMPs 
should. Moreover, the results of analysis of the outcome of a direct detection experiment would depend on the value 
of the escape speed chosen, which is quite uncertain (see, e.g., Ref. [HI). 
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FIG. 6: 90% C.L. allowed region in the WIMP mass vs. WIMP-nucleon spin-independent cross section plane implied by the 
DAMA modulation signal, for our best-fit halo model with parameters rt = 120 kpc, a = 300 km s~^ and Pdm,q ~ 0.2 GeV/ cm^ . 
Also shown are 90% C.L. upper limits on cross section as a function of WIMP mass as implied by the null results of other 
experiments, as indicated, again for the same King model halo parameters. 




(GeV) (GeV) 

FIG. 7: 90% C.L. allowed region in the WIMP mass vs. WIMP-neutron (op = 0) (left panel) and WIMP-proton (a„ = 0) (right 
panel) spin-dependent cross section plane implied by the DAMA modulation signal, for our self-consistent truncated isothermal 
(King) model of Milky Way's DM halo with parameters rt = 120 kpc, a = 300 kms^^ and Pdm,© ~ 0.2 GeV/ cm'^. Also shown 
are 90% C.L. upper limits on the cross sections as a function of WIMP mass as implied by the null results of other experiments, 
again for the same King model halo parameters. 



In this paper we have presented a self-consistent model of the finite-size dark halo of the Galaxy with its phase 
space distribution function described by a truncated isothermal model ("King" model). We have also included the 
gravitational influence of the observed visible matter on the structure of the DM halo in a self-consistent manner. In 
this model the velocity distribution function of the WIMPs constituting the halo is non-Maxwcllian, with a cut-off at a 
maximum velocity that is self-consistently determined by the model itself. We have determined the parameters of our 
halo model by a fit to a recently determined circular rotation curve of the Galaxy that extends up to a Galactocentric 
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distance of ^ 60 kpc [29[ . A noticeable feature of this rotation curve is that it dechnes with the Galactocentric radius 
beyond ~ 10 kpc. This imposes strong constraints on the (truncation) radius, density and mass of the DM halo 
and hence on the total mass of the Galaxy. Specifically, for values of local DM density in the range Pdm.o = 0.2 
- 0.4GeV/cm'^, reasonable fits to the rotation curve data require the truncation radius of the halo to be in the 
corresponding range rt ~ 120 - 80 kpc, in all cases restricting the total mass of the Galaxy (including its DM halo) to 
relatively low values of Af Galaxy ^ 3 x 10^^ Mq. Note, however, that the rotation curve estimated in Ref. [2^, which 
is used in determining the parameters of our halo model, was derived assuming a value of w 220 km s"-'^ for the local 
standard of rest. On the other hand, as suggested in Ref. [2^, the value of the local standard of rest may have to be 
revised upward to a value of Vc.q ~ 250 km s^^. This would require appropriate scaling up of the rotation curve data 
of Ref. [29|, which in turn will require an appropriately higher value of Pdm,© a-nd a correspondingly higher value of 
the total mass of the Galaxy. 



Interestingly, we find that the upper limits on the relevant WIMP-nuclcon interaction cross section implied by the null 
results of the direct detection experiments are primarily determined by the chosen value of the local DM density Pdm,0 
— scaling roughly inversely with the value of Pdm,© — and are relatively less sensitive to the other parameters of the 
model such as the truncation radius and total mass of the halo. For our best-fit self-consistent lowered isothermal 
halo model with parameter values Pdm,o = 0.2 GeV/ cm^ and rt = 120 kpc, that provides the best fit to the rotation 



curve data of Ref. [29[, the null result of the CDMS-II experiment [ll|, for example, gives a 90% C.L. upper limit on 
the WIMP-nucleon spin-independent (SI) interaction cross section of ~ 5.3 x 10~^pb at a WIMP mass of 71 GcV. 



Concerning the issue of the compatability of the claimed positive signal reported by the DAMA collaboration with the 
null results from the other experiments, we find, within the context of a simple analysis procedure using the original 
2-bin DAMA annual modulation data, that there exist regions of the WIMP mass vs. WIMP-nucleon cross section 
with small WIMP masses typically in the range '--^ 2 - 16 GeV, within which the DAMA's claimed annual modulation 
signal is consistent with the null results of other experiments. While this strengthens the possibility of a low mass 
WIMP as a DM candidate — a possibility also indicated by several earlier analyses done using the SHM — a more 
rigorous analysis using DAMA bins over smaller intervals must be performed before a more definitive conclusion can 
be reached in this regard. Very recently, the CoGeNT collaboration [s^l has also reported an excess of low (nuclear 
recoil) energy events, which may also be pointing towards a low-mass WIMP DM candidate. We shall discuss these 
issues within the context of our self-consistent model of the Galaxy's dark halo in a work now in progress. 
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